Probing nuclear symmetry energy with the sub-threshold pion production 
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Within the framework of semiclassical Boltzmann-Uehling-Uhlenbeck (BUU) transport model, 
we investigated the effects of symmetry energy on the sub-threshold pion using the isospin MDI 
interaction with the stiff and soft symmetry energies in the central collision of 48 Ca + 48 Ca at the 
incident beam energies of 100, 150, 200, 250 and 300 MeV/nucleon, respectively. We find that 
the ratio of n~ /-k + of sub-threshold charged pion production is greatly sensitive to the symmetry 
energy, particularly around 100 MeV/nucleon energies. Large sensitivity of sub-threshold charged 
pion production to nuclear symmetry energy may reduce uncertainties of probing nuclear symmetry 
energy via heavy-ion collision. 
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The density dependence of the nuclear symmetry en- 
ergy is not only important for nuclear physics, but also 
crucial for many astrophysical processes, such as the 
structure of neutron stars and the dynamical evolution 
of proto-neutron stars [l| . Though considerable progress 
has been made recently in determining the density de- 
pendence of the nuclear symmetry energy around the 
normal nuclear matter density from studying the isospin 
diffusion in heavy-ion reactions at intermediate energies 
@~Bl> much more work is still needed to probe the high- 
density behavior of the nuclear symmetry energy. For- 
tunately, heavy-ion reactions, especially those induced 
by radioactive beams, provide a unique opportunity to 
constrain the EOS of asymmetric nuclear matter, and 
a number of such observables have been already identi- 
fied in heavy-ion collisions induced by neutron-rich nu- 
clei, such as the free neutron/proton ratio [f|, the isospin 
fractionation @, @], the neutron-proton transverse dif- 
ferential flow [8| , the neutron-proton correlation function 
0,t/3He[lo|,the isospin diffusion [ll[ , the proton differ- 
ential elliptic flow [12|. Currently, to pin down the sym- 
metry energy, the National Superconducting Cyclotron 
Laboratory (NSCL) at Michigan State University, the 
Geselschaft fuer Schwerionenforschung (GSI) at Darm- 
stadt, the Rikagaku Kenkyusho (RIKEN, The Institute 
of Physical and Chemical Research) of Japan, as well as 
the Cooler Storage Ring (CSR) in Lanzhou, are planning 
to do related experiments. Some of related information 
can be find via Ref. jl3| . 

Recently, pion production in heavy-ion collisions has 
attracted much attention in the nuclear physics commu- 
nity 0, [H, [Til [2l|. One important reason for this is 
that pion production is connected with the high-density 
behavior of nuclear symmetry energy, especially around 
pion production threshold. The latter is crucial for 
understanding many interesting issues in both nuclear 
physics and astrophysics. And several hadronic trans- 
port models have quantitatively shown that ir~ /ir + ra- 
tio is indeed sensitive to the symmetry energy, especially 
around pion production threshold [19j and above pion 



production threshold [24|. In the previous studies, pion 
production mainly from energetic heavy-ion collisions, in 
which pions are mainly produced above their threshold 
energy. What is the case of effects of nuclear symmetry 
energy on charged pion ratio via sub-threshold produc- 
tion? To answer this question, we studies the effects of 
nuclear symmetry energy on sub-threshold pion produc- 
tion. We select 48 Ca+ 48 Ca as the reaction system due 
to it large asymmetry. Based on the isospin-dependent 
Boltzmann-Uehling- Uhlenbeck (IBUU) transport model, 
we studied the effects of the symmetry energy in the cen- 
tral reaction 48 Ca+ 48 Ca and find that the ratio of ir~ /ir + 
of sub-threshold charged pion production, compared with 
above threshold case, is particularly sensitive to the sym- 
metry energy. 

Our present work is based on the semi-classical trans- 
port model IBUU04, in which nucleons, A and N* res- 
onances as well as pions and their isospin-dependent 
dynamics are included. We use the isospin-dependent 
in-medium nucleon-nucleon (NN) elastic cross sections 
from the scaling model according to nucleon effective 
masses. For the inelastic cross sections we use the ex- 
perimental data from free space NN collisions since at 
higher incident beam energies the NN cross sections have 
no evident effects on the slope of neutron-proton differen- 
tial flow. The total and differential cross sections for all 
other particles are taken either from experimental data 
or obtained by using the detailed balance formula. The 
isospin dependent phase-space distribution functions of 
the particles involved are solved by using the test-particle 
method numerically. The isospin-dependence of Pauli 
blockings for fermions is also considered. The isospin- 
dependence of Pauli blockings for fermions is also con- 
sidered. More details can be found in Refs. [3, IM Il6l - fl9| . 
In the present studies, the momentum-dependent single 
nucleon potential (MDI) adopted here is: 
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In the above equation, S — {p n — p P )/(p n + Pp) is the 
isospin asymmetry parameter, p = p n + p p is the baryon 
density and p n , p p are the neutron and proton densi- 
ties, respectively, r = l/2(— 1/2) for neutron (pro- 
ton) and r ^ r', a = 4/3, / T (r,p) is the phase- 
space distribution function at coordinate r and momen- 
tum p. The parameters A u (x), Ai(x), B,C TiT , C TiT > and 
A were set by reproducing the momentum-dependent 
potential U(p, 5, p, r) predicted by the Gogny Hartree- 
Fock and/or the Brueckner-Hartree-Fock calculations. 
The momentum-dependence of the symmetry potential 
steams from the different interaction strength parame- 
ters C TtT i and C TtT for a nucleon of isospin r interacting, 
respectively, with unlike and like nucleons in the back- 
ground fields, more specifically, C un iike = —103.4 MeV 
while Cuke = —11.7 MeV. The parameters A u (x)a.nd 
Ai{x) depend on the x parameter according to Au{x) = 
-95.98- ir -2#r an&Adx) = -120.57+z- 2 ^-. The satura- 
tion properties of symmetric nuclear matter and the sym- 
metry energy of about 32 MeV at normal nuclear matter 
density po — 0.16 fm -3 . The incompressibility of sym- 
metric nuclear matter at normal density is set to be 211 
MeV. According to essentially all microscopic model cal- 
culations, the EOS for isospin asymmetric nuclear matter 
can be expressed as 



E(p, 6) = E(p, 0) + E sym (p)S 2 + 0{6 4 ), 



(2) 



where E(p, 0) and E sym (p) are the energy per nucleon of 
symmetric nuclear matter and nuclear symmetry energy, 
respectively. For a given value x, with the single parti- 
cle potential U(p, S, p,r), one can readily calculate the 
symmetry energy E sym (p) as a function of density. 

The main reaction channels related to pion production 
and absorption are 
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where R denotes A or N* resonances. In the present 
work, we use the isospin-dependent in-medium reduced 
NN elastic scattering cross section from the scaling 
model according to nucleon effective mass [ll| to study 
the effect of symmetry energy on pion production. As- 
suming in-medium NN scattering transition matrix is 
the same as that in vacuum, the elastic TV TV scattering 
cross section in medium <r^^ mm is reduced compared 
with their free-space value ujj^ e by a factor of [2(| 
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where and p* NN are the reduced masses of the col- 
liding nucleon pair in free space and medium, respec- 
tively. For in-medium NN inelastic scattering cross sec- 
tion, even assuming in-medium NN — > NR scattering 
transition matrix is the same as that in vacuum, the den- 
sity of final states Dj, of NR is very hard to calculate 
due to the fact that the resonance's potential in matter 
is presently unknown. The in-medium NN inelastic scat- 
tering cross section is thus quite controversial. Because 
the purpose of present work is just study the effect of 
symmetry energy on pion production and charged pion 
ratio, to simplify the question, for the A^iV inelastic scat- 
tering cross section we use the free A^iV inelastic scatter- 
ing cross section. The effective mass of nucleon in isospin 
asymmetric nuclear matter is 
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From the definition and Eq. (p}, we can see that the ef- 
fective mass depends not only on density and asymmetry 
of medium but also the momentum of nucleon. 
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FIG. 1: (Color online) Density dependence of nuclear sym- 
metry energy with parameters x — 1,-1, respectively. 

Fig. Q] shows the density dependence of nuclear sym- 
metry energy with parameter x = 1,-1, respectively. 
As discussed in the previous part, the single particle 
used has an x parameter, different specific x parame- 
ter denotes different density dependent symmetry energy. 
For the central reaction 48 Ca+ 48 Ca, the maximal den- 
sity reached is about 1.5 ~ 2 times saturation density as 
shown in Fig. [2] From Fig. [1] we can also see that the 
low density behaviors of nuclear symmetry energy sepa- 
rate from each other with different x parameters. At the 
saturation point there is a cross and then they separate 
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FIG. 2: (Color online) Maximal baryon density reached in 
the central reaction 48 Ca+ 48 Ca with different incident beam 
energies. 



from each other again. At lower densities, the value of 
symmetry energy of x — is lower than that of x = 1. 
But at high densities, the value of symmetry energy of 
x = is higher than that of a; = 1. From Fig. [3J we 
can see that the maximal baryon density reached in the 
central reaction 48 Ca+ 48 Ca increases with the incident 
beam energy. At 100 MeV/nucleon, the maximal baryon 
density reaches about 1.5 times saturation density and at 
300 MeV/nucleon, the maximal baryon density reaches 
about 2 times saturation density. Therefore the ratio 
of 7r~ /ir + of sub-threshold charged pion production still 
mainly reflects the high density behavior of nuclear sym- 
metry energy. 

Fig. [3] shows the kinetic energy distribution of the 
7r~/7r + ratio using the MDI interaction with x = 1 and 
-1 for the central collision of 48 Ca + 48 Ca at the incident 
beam energies of 150 MeV/nucleon. It is seen that at pion 
kinetic energies of 30 ~ 120 MeV, the ir~ /ir + ratio is very 
sensitive to the symmetry energy. The soft symmetry 
energy (x= 1) gives larger tt~ /ir + ratio whereas the stiff 
symmetry energy (x = —1) gives smaller tt~Jjt + ratio. 
This is consistent with the previous studies [li], Oil IHJ 
for the neutron- rich reactions. And we can see that the 
sensitivity of ir~ /ir + ratio to the high density behavior 
of symmetry energy is quite large, about 80%, is quite 
sensitive to cases of pion production at higher incident 
beam energies [22| • The sub-threshold charged pion pro- 
duction is thus a very sensitive probe of the symmetry 
energy. Fig. 2] is the case at 300 MeV/nucleon incident 
beam energy. Again we clearly see that at pion kinetic 
energies of 30 ~ 120 MeV, the tt~ /tt + ratio is sensitive to 
the symmetry energy. The soft symmetry energy (x = 1) 
gives larger tt~ /tt + ratio whereas the stiff symmetry en- 
ergy (x — —1) gives smaller tt~ /n + ratio. Sensitivity of 



FIG. 3: (Color online) Kinetic energy distribution of the 
7r~ /tx + ratio using the MDI interaction with x = 1 and -1 
for the central collision of 48 Ca + 48 Ca at the incident beam 
energies of 150 MeV/nucleon. 
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FIG. 4: (Color online) Kinetic energy distribution of the 
7r~ /tx + ratio using the MDI interaction with x = 1 and -1 
for the central collision of 48 Ca + 48 Ca at the incident beam 
energies of 300 MeV/nucleon. 



charged pion ratio to the symmetry energy at higher in- 
cident beam energy 300 MeV/nucleon is smaller than the 
case at lower incident beam energy 150 MeV/nucleon. In 
the heavy-ion collisions at intermediate energies, gener- 
ally speaking, mean-field effect and collision effect com- 
pete each other. The effect of mean-field increases with 
the decreasement of incident beam energy. Thus at lower 
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incident beam energies effects of the symmetry energy 
(which is the isovector part of nuclear mean-field poten- 
tial) are larger than that with higher incident beam en- 
ergies. From Fig. 2J we can also see that there is an evi- 
dent Coulomb peak [l9j], and around the Coulomb peak 
7r~/7r + ratio is more sensitive to the symmetry energy. 
At lower incident beam energy 150 MeV/nucleon, how- 
ever, there is no Coulomb peak at at pion kinetic energies 
of 30 ~ 120 MeV. 
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FIG. 5: (Color online) Excitation function of the tt~ /n + ratio 
using the MDI interaction with x — 1 and -1 for the central 
collision of 48 Ca + 48 Ca at the incident beam energies of 100, 
150, 200, 250 and 300 MeV/nucleon, respectively. 

Fig. [5] shows the excitation function of the ir~ /ir + ra- 
tio using the MDI interaction with x = 1 and -1 for the 
central collision of 48 Ca + 48 Ca at the incident beam 
energies of 100, 150, 200, 250 and 300 MeV/nucleon, re- 
spectively. From Fig. [5] we can clearly see that sensitiv- 



ity of charged pion ratio tt~ /ir + reaches the maximum at 
the lower incident beam energy 100 MeV/nucleon, about 
100% sensitive to the symmetry energy. The trend of the 
sensitivity of the charged pion ratio 7r~ /ir + to symmetry 
energy decreases with the incident beam energy is consis- 
tent with the previous studies [l4| • Note here that at the 
incident beam energy region studied here, the soft sym- 
metry energy always corresponds large value of 7r~ /tt + 
ratio and the stiff symmetry energy corresponds relative 
small value of n~ /n + ratio (23|. However, in Ref. fzij . 
Feng et al. systematically investigated the pion produc- 
tion in heavy-ion collisions in the region of below 1 ^4GeV 
energies by using the ImlQMD model. They found that 
the excitation functions of the tt~ /tt + ratio increases with 
the stiffness of the symmetry energy. This is inconsistent 
with our studies qualitatively and the reasons are needed 
to be clarified. 

In summary, based on the semiclassical Boltzmann- 
Uehling-Uhlenbeck (BUU) transport model, we studied 
the effects of symmetry energy on the sub-threshold pion 
production in the central reaction 48 Ca+ 48 Ca. We find 
that the ratio of ir~ of sub-threshold charged pion 
production is particularly sensitive to the symmetry en- 
ergy with decrease of incident beam energy of heavy-ion 
collision. The highly sensitive charged sub-threshold pion 
ratio to the symmetry energy may reduce the uncertain- 
ties of probing nuclear symmetry energy via heavy-ion 
collision. 
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